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Abstract 

In this paper we present the complete two-loop vertex corrections to scalar 
and pseudo-scalar Higgs boson production for general colour factors for the 
gauge group SU(N) in the limit where the top quark mass gets infinite. We 
derive a general formula for the vertex correction which holds for conserved 
and non conserved operators. For the conserved operator we take the elec- 
tromagnetic vertex correction as an example whereas for the non conserved 
operators we take the two vertex corrections above. Our observations for 
the structure of the pole terms and in two loop order are the 

same as made earlier in the literature for electromagnetism. However we also 
elucidate the origin of the second order single pole term which is equal to 
the second order singular part of the anomalous dimension plus a universal 
function which is the same for the quark and the gluon. 
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1 Introduction 



One of the crucial tests of the standard model will be the discovery of the 
Higgs boson at the LHC or at the TEVATRON. Its discovery or its absence 
will shed light on the mechanism how particles acquire their mass. It will also 
answer questions about supersymmetric extensions of the standard model or 
about compositeness of the existing particles including the Higgs boson. In 
this paper we will concentrate on Higgs production where the lowest order 
reaction proceeds via the gluon-gluon fusion mechanism. In the standard 
model the scalar Higgs boson H couples to the gluons via heavy quarks 
among which the top-quark is the most prominent one. This also holds for 
the pseudo-scalar Higgs boson A provided tan /3 is small where (3 denotes the 
mixing angle in the Two-Higgs-Doublet model (2HDM). The basic produc- 
tion mechanism is that two incoming gluons are coupled to the Higgs boson 
through a triangular top-quark loop [1]. The first order QCD corrections to 
this process [2] are already very complicated in particular the virtual cor- 
rections which involve massive two-loop integrals. However fortunately one 
discovered that one could make a simplification for the total cross section 
[2], [3]. The latter is also possible for differential distributions provided the 
transverse momentum is not too large [px ^ '^itt-h) [4], [5], [6], [7]. This 
is achieved by taking the large top-quark mass limit rrit ^ oo. In this case 
the Feynman graphs are derived from an effective Lagrangian describing the 
direct couphng of the (pseudo-) scalar Higgs boson to the gluons. Therefore 
the triangular top-quark loop is replaced by a simple vertex and the treat- 
ment of the total cross section is the same as for the Drell-Yan process. In 
this paper we concentrate on the one and two-loop vertex QCD corrections 
to the total cross section for Higgs boson production. For the two-loop elec- 
troweak corrections we refer to [8] and [9]. Although the two-loop virtual 
contributions for scalar Higgs production in the heavy top-quark limit were 
presented in [10], the result was not decomposed in the various colour fac- 
tors . This is important for the general structure of the vertex function as 
will be pointed out later on. An explicit expression for the two-loop virtual 
corrections to the pseudo-scalar Higgs boson is not given anywhere although 
it is indirectly contained in [11], [12]. In this paper we explicitly present the 
two-loop virtual contributions to scalar as well as pseudo-scalar Higgs boson 
production for general colour factors in SU(N). 

The paper will be organized as follows. In section 2 we give an outhne 
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of the renormalization constants which are involved in the calculation of the 
two-loop virtual corrections to scalar and pseudo-scalar Higgs production. 
In section 3 we present the calculation and compute the contribution to the 
flavour dependent part of the coefficient function in the total cross section 
for (pseudo-) scalar Higgs production. In section 4 we discuss the structure 
of the vertex correction w.r.t. the pole terms and present the finite terms. In 
particular we find that also the single pole term can be predicted so that the 
vertex and other two-loop radiative corrections can be expressed in a unique 
form. For completeness we also compare it with the electromagnetic form 
factor present e.g. in the Drell-Yan process. 
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2 Application of the effective Lagrangian ap- 
proach to Higgs production. 

In the large top-quark mass limit the Feynman rules (see e.g. [13]) for scalar 
Higgs production (H) can be derived from the following effective Lagrangian 
density 



£H ^ = Gh <^^{x) 0{x) with 0{x) = 
whereas pseudo-scalar Higgs (A) production is obtained from 



-Gl^{x)G'^'^'^{x), (2.1) 



GaOi(x) + Ga02(x] 



with 



Oi(x) = -^e^.A.G'rG'f(x), 



^ i=i 



(2.2) 



where $^(x) and $^(x) represent the scalar and pseudo-scalar fields respec- 
tively and Uf denotes the number of light flavours. Furthermore the gluon 
field strength is given by G^" and the quark field is denoted by g^. The fac- 
tors multiplying the operators are chosen in such a way that the vertices are 
normalised to the effective coupling constants Gh_, Ga and Ga- The latter 
are determined by the top-quark triangular loop graph, including all QCD 
corrections, taken in the limit — > oo which describes the decay process 
B ^ g + g with B = H, A namely 

Gb = -2^/'a,(/x^)Gf rBFB(rB)CB(^a,(/x^),^j , 



Ga , 



and as{n1) ^he renormalized coupling constant defined by 



47r 



(2.3) 



(2.4) 
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where as(yU^) is the running couphng constant and fir denotes the renormal- 
ization scale. Further Gp represents the Fermi constant and the functions 
Fb are given by 

Fnir) = 1 + (1 - r) /(r) , FA(r) = /(r) cot /? , 
4 m? 



fir) 



arcsm 



for r > 1 , 



/w = -J 



In 



1 + ^ 



TT 2 



for T < 1 , 



(2.5) 



where cot f3 denotes the mixing angle in the Two-Higgs-Doublet Model. Fur- 
ther m and rrit denote the masses of the (pseudo-) scalar Higgs boson and 
the top quark respectively. In the large mrlimit we have 



lim Fnir) = — , 

r— »oo o T 



lim Fa(t) = - cot/?. 

r— »oo T 



(2.6) 



The coefficient functions Cb originate from the corrections to the top-quark 
triangular graph provided one takes the limit — > oo. We have presented 
the Born level couphngs Gb in Eq. (2.3) for general rut for on-shell gluons 
only in order to keep some part of the top-quark mass dependence. This is 
an approximation because the gluons which couple to the H and A bosons 
via the top-quark loop in partonic subprocesses are very often virtual. The 
virtual-gluon momentum dependence is neither described by Fb{t) nor by Cb- 
However for total cross sections the main contribution comes from the region 
where the gluons are almost on-shell so that this approximation is better 
than it is for differential cross sections with large transverse momentum. 
The coefficient functions are computed up to order in [14], [15] for the H 
and in [16] for the A. They read as follows 



V rntj 



5Ga — 3 Ci 



2 ^ 



CATf + {lC^ 



5 



mj \ 9 mi 



(2.7) 



CA(a,(/.^),^j =1, (2.8) 

where af"^ is presented in the five-flavour number scheme. Notice that the 
coeflicient function in Eq. (2.7) is derived for general colour factors of the 
group SU(N) from Eq. (6) in [15]. These factors are given by 

A^^ - 1 1 

Ca = N , Cf= , Tf = -. (2.9) 

A ' 2A^ 2 ^ ^ 

Notice that Tf is also incorporated into in Eq. (2.3) where it is set to the 
value Tf = 1/2. Using the effective Lagrangian approach we will calculate 
the total cross section of the reaction 

//i(Pi) + 7/2(^2) ^ B(g) +' X' , (2.10) 

where Hi and H2 denote the incoming hadrons and X represents an inclusive 
hadronic state. The total cross section is given by 

° ^1 a,b=q,q,g ''^ ''^l^^ 

(X Tfl^ \ 
,— , B = H,A, 
Xi X2 jJ.^ I 



with x^-^ , S^{Pi + P2f , q^^m\ (2.11) 

where the factor 1/{N^ — 1) originates from the colour average in the case of 
the local gauge group SU(N). Further we have assumed that the (pseudo-) 
scalar Higgs boson is mainly produced on-shell i.e. ~ m^. The parton 
densities denoted by faiVyfJ^"^) {ci,b = q,q,g) depend on the mass factoriza- 
tion/renormalization scale /x. The same scale also enters the coefficient func- 
tions Aab,B which are derived from the partonic cross sections {z — w? / s) 

^,^1 = -y j^,S^iq^-m^)Tab,B{(l^Pi,P2), (2.12) 



where the incoming partons a and h carry momenta pi and p2 respectively. 
They are related to the hadron momenta by 



S={P1+P2) , S = XiX2S, . 

s 

The amplitude Tab,B can be written as 

TabAQ,Pi,P2) = Glld'ye"^y{a,b\0(y)Oma,b), 

TabAq:Pl:P2) = / A ' M^, ^'1 (^A d (y) + O2 (y)) 



(2.13) 
(2.14) 



X (Ga Oi(0) + Ga 02(0)) \a, b) . (2.15) 

The expressions above for the amplitude T^b are similar to those given for 
the Drell-Yan process except that the conserved electroweak currents are 
replaced by the operators O and Oi,02. The latter are not conserved so 
that they acquire additional renormalization constants defined by 

0{y) = Zod{y) , 0,{y) = Z,,dj{y). (2.16) 

where the hat indicates that the quantities under consideration are unrenor- 
mahzed. Insertion of the above equations into Eqs. (2. 14), (2. 15) leads to the 
renormalized expressions 

TabAQ,PuP2) = GlZl j Ae^^-M«,&|O(2/)O(0)|a,6), (2.17) 



^11 + ^21 + 2 Ga Ga Zii Z- 



A '^A ^11 ^21 



X {a, b\0,{y) Oi(0)|a, b) + {Gl Zf^ + Gl Zl^ + 2 Ga Ga ^12 ^22 



X (a, b\62{y) C>2(0)|a, 6) + G^ ^12 + G\ Z22 ^21 + Ga Ga 



X (Zii Z22 + ^12 ^21) [ (a, 6|Oi(y) 62(0) + ^(y) Oi(0)|a, 6) 



(2.18) 



Since we are interested in the one- and two-loop correction to the subprocess 
B ^ g + g the above formula simplifies and it becomes 

W?>Pi>P2) = GlZ'o ld'ye'^y{g,g\d{y)d{0)\g,g), (2.19) 



Gl Zl, {g, g\0^{y) 0,(0)19, 9) + {Gl ^12 



+Ga Ga {g, g\d,{y) 02(0) + O^iy) dM\9, 9) 



.(2.20) 



The operator renormalization constants depend on the regularization scheme 
in particular on the prescription for the 75-matrix and the Levi-Civita tensor 
in Eq. (2.2). The computation of T,,(;.h will be carried out by choosing n- 
dimensional regularization where in the case of Tgg^^ we adopt the HVBM 
prescription [17], [18] for the 75-matrix. For this choice the contraction of 
the Levi-Civita tensors proceeds as 



^/ili^lAicri ^ 



^^2 §\2 ^<T2 

"^1 "^^i "^/Xi 

Kl 511 

K ^A? 

<^ 5^1 K\ 



Ml 



(2.21) 



where all Lorentz indices are taken to be n-dimensional. To facilitate the 
calculation one can replace 75 in Eq. (2.2) by (see [19]) 



75 



24 



(2.22) 



which is equivalent to the HVBM scheme. Choosing the MS subtraction 
scheme the renormalization constant corresponding to the operator O be- 
comes [20] 



£ 



2 

£ 



+ 



(2.23) 



where denotes the spherical factor characteristic of n-dimensional regular- 
ization. It is defined by 



= exp <{ - 



7e 



ln47r 



}■ 



(2.24) 



8 



The lowest order coefficients (3o and /3i originate from the beta-function given 
by 



11 „ 4 „ „ 34 



20 



(2.25) 

The operator renormahzation constants corresponding to Oi and O2 are com- 
puted in [21] and they read 



Zii = = 1 + asii^i) Se-Po + aM) S; 

£ 



+ 



(2.26) 



where Z^^ denotes the couphng constant renormahzation factor defined by 

as = Z^^ a^ifil) . (2.27) 
The remaining constants are 

Z21 = 0, (2.28) 
1 



Zu = as{ijL^)S^- -6Cf 

£ L 

^22 — ^MS ^5 > 



(2.29) 

(2.30) 



where Z^g and Z| are the constants characteristic of the HVBM scheme. 



They are given by 



—^A^F — J- S^F — 



Zl = l + a,(/x^) 



4Cf 



+ a,(/ij 



22 



31 



F ^ C'a + n/ T/ CiT' 



(2.31) 



The latter renormahzation constant is determined in such a way that the 
Adler-Bell-Jackiw anomaly [22] 

O^iy) = ^as{[il)n^T^O^{y) , (2.32) 
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is preserved in all orders in perturbation theory according to the Adler- 
Bardeen theorem [23]. 
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3 Two-loop vertex correction to the process 
g + g^H,A. 

In this section we present the calculation of the (pseudo-) scalar Higgs vertex 
corrected up to two loops. The scalar Higgs vertex was presented in [10] 
but no separation in colour factors was made. An explicit formula for the 
pseudo-scalar vertex was not shown but it indirectly appears in the coefficient 
functions of Higgs boson production in [11] and [12]. 

The one loop vertex correction is very easy. In the case of Higgs produc- 
tion the characteristic graphs are shown in Fig. 1. For the operators 0{y) 
and Oi{y) all the solid lines represent gluons and there are only two types 
of graphs. The two-loop vertex correction, which is less trivial, can be cal- 
culated in various ways. The characteristic graphs are shown in Fig. 2. The 
solid lines represent gluons, quarks as well as ghosts (Feynman gauge). Using 
the program FORM [24] each graph can be reduced into scalar integrals of 
the type 

^ r dr%_ r dl'k2 {k ■ piY {k^ ■ . ^ 

^n3.-3. J ^^^y] (2^)„ ''^ ' 

Da = (/C2 - Pi? , D, = {k2 - P2f , De = kj , 

DT^{ki-pif, Ds^{ki-p2f. (3.1) 

The scalar integrals in the case of the Higgs boson vertex are ordered as 
follows. For six propagators we have i+j — A. For five propagators i+j — 3. 
For four propagators i+j = 2. For three propagators i+j — 1. The two and 
one propagator integrals are all zero. The calculation of these scalar integrals 
can be performed in various ways. For the electromagnetic vertex correction 
in QCD it has been done in [25], [26], [27]. In [25] one has chosen a Feynman 
parameter set in a judicious way so that after each integration one is still left 
with a linear integral over the next parameter. In [26] one used integration by 
parts according to [28] . This trick eliminates one propagator from the integral 
in Eq. (3.1). However it does not work for the second Feynman graph in 
Fig. 2 (the crossed box graph). In [27] the scalar integrals were calculated 
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Figure 1: General structure of the one- loop vertex correction g-g-B {B — 
H,A). 



using dispersion techniques developed in [29] . In the latter reference one has 
used the Cutkosky rules [30] to compute the scalar integrals. Here one cuts 
the graph in all possible ways to obtain the imaginary part. The real part 
is obtained via a dispersion relation. Fortunately for the electromagnetic 
vertex correction i + j < 3 in Eq. (3.1) but for the vertex correction to 0{y) 
and Oi{y) in Eq. (2.2) one also gets contributions for i + j = 4. In [10] 
one employed an algorithm in [31] which relates 1-loop integrals with n + 1 
external legs to I + 1-loop integrals with n external legs. In the case under 
consideration massless two-loop vertex functions are mapped onto massless 
three-loop two-point functions. We however adopt the method in [29] and 
calculate the imaginary part of all graphs using the Cutkosky rules. In the 
case that all particles are massless, except for the external Higgs boson leg, 
the dispersion integrals are trivial. Since all graphs, except for the second 
graph in Fig. 2, are analytic we have made a computer program in FORM 
[24] which provides us with all scalar integrals of the type in (3.1). For more 
details about the real and imaginary parts of the scalar integrals we refer to 
[29]. The result for the scalar Higgs vertex is obtained from Eq. (2.19) and 
it reads up to order 

Fniq') (27r)^5(g + pi -p2) = / Ae^^'^ (^(P2)|0(y)|^(pi)) , (3.2) 
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where Fniq^) for < is equal to {n — A-\-e) 



2 \ e/2 



J + C(2) + s(l-k(3)) 



+ a,S, i — 



32 44 
+ 



/134 , ,x 1 / 136 ^, , 50 ^, ,x 1 
-(^ + 4 C(2)) + ( - 11 C(2) + - C(3)) - 



+ 



9 V £2 ■ V 27 

5861 67 ^,^, 11,,,, 21 

\ C 2 H C(3 C 2 

162 6 ^ 9 ^ 5 ^ ^ ^ 



14 




1 1616 



81 



fc(2)-^C(3)- 



[ 4 67 



(3.3) 



where hg (see Eq. (2.4)) is the bare couphng. The result above agrees with 
Eqs. (16) and (17) in [10] provided we put Ca = 3 and Cp = 4/3 (Notice 
that in [10] n = 4 — 2 e). The pseudo-scalar Higgs vertex can be derived from 
Eq. (2.20). 



FA{q^){2TiY5{q + p,-p2) 



d'ye'^-y\Zn (<7(P2)|Oi(y)|<7(pi)) 



{9ip2)\02{y)\gipi))\, 



(3.4) 



The matrix element {g{p2)\02{y)\g{pi)) = C ag but Z12 ~ asB~^ (see Eq. 
(2.29)). This contribution has to be added because we use the HVBM scheme 
for the 75-matrix. The extra term equals (see Fig. 3) 



24 



-36 



(3.5) 
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In the equation above we have multiphed by an extra factor Ss {—(f/jJ-'^Y^'^ 
to restore the anti-commutativity of the 75-matrix. In a scheme where the 
75-matrix anti-commutes with the other gamma matrices this extra factor 
would be automatically there. This term has to be added to F^iq^) and the 
result up to order becomes {(f < 0) 



4 + C(2) + £(-6-^C(3 



/7_1^(2) + -C'(2 
V 2 ^ 80^ ^ 



44 
"3^ 



2 

/422 
v~9~ 




27 
21 



3835 115 , 11 ^ .2/ X 
+^ + — C(2) + - C(3) - - C'(2) 



81 



6 



9 



32 



fc(3))i 



16 



40 



3£3 



+ 



/212 



+ 4C(2))1 



3182 



-^C(2)-^C(3)- 



+nfTfCF 



142 



+ 16C(3) 



(3.6) 



This result is new. It only indirectly appears in the 6{1 — z) part of the total 
cross section for pseudo-scalar Higgs boson production in [11] and [12]. We 
can now compute the part due to virtual and soft gluons which is proportional 
to S{1 — z) for general colour factors as in [32]. Continuing the above form 
factors to > and adding the soft gluon cross section we can write the 
soft plus virtual gluon contribution to the coefficient function in Eq. (2.11) 
as 



aS+V 



'gg,U > gg,B 

(3.7) 

Starting with scalar Higgs boson production we omit the term which is pro- 
portional to as it is available in [32] . The remaining order part of Eq. 



16 



(3.7) equals 



. (2),S+V _ ^ rj. ^ 



16 



Vo{z) In^ 



m 



64 _ , , 160 , 



+5{1 - z) (l6 + - C(2)) I In i^—j +-V2{z)-— V, 
/448 64 ,,^,x ^ , ^ , / 160 160 



-tC(3)) 



+71/7/(7^5(1-2) 



, /m^\ 134 



(3.8) 



The result above is in agreement with [11] and [33] for Ca = 3 and Cp = 4/3. 
For the pseudo scalar case we take the difference with the scalar Higgs boson 
coefficient function. Moreover we add the last term in Eq. (2.20) which 
equals 

I d^y e'^^-y Ga Ga {g,g\d,{y) 62(0) + 62(1/) Oi(0) \g, g) . (3.9) 

This term is proportional to O12 in the formula below and it originates from 
Eq. (2.3). Ignoring the term proportional to C\ which is given in [32] the 
remaining order term equals 



A(2).S+V 



Ca TfUfSil — z 
+CFTfnf5{l-z) 



- - In 
3 



m 



- 8 In 



m 



11^ 



-50 



+O12 24 In 



- 12 



(3.10) 



which agrees with [11] and [12] for Ca — 3 and Cp = 4/3. 
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4 General structure of the vertex correction 



In this section wc will study the general structure of the vertex correction as 
predicted in [34] and [35]. In these two papers the behaviour of the vertex 
corrections w.r.t. ln(— was derived. Following Eq. (6.21) in [35] the 
behaviour is as follows 



dhi{-q^/jj?) 



^ lK,^^(as{^^l)) + ^ Gii(as{q^)) 



1 



+Kii ( -,as(/x 



where 



and 



,3/',,2n ,(2) 



(4.1) 



(4.2) 



Gii 



■ (0) 

Tii - zi 



(1) 

Tii -Z2 



(2) 



^3 



+ 



The quantities ^K,ii and Ga {i — q, g) occur in the diagonal kernel 

1 



1 + 



lK,v 



1 



5(1 -x) 



(4.3) 



(4.4) 



and they represent those parts of the splitting functions which are propor- 
tional to the distributions (1/(1 — a;))+ and 5(1 — x) respectively provided 
they are expressed in the MS-scheme. They have been calculated up to three 
loop order in [36] . The constants zi occur when the operator is not conserved 
which happens in the case of Higgs boson vertex corrections. The constants 
show up in the operator renormalization constants Zq in Eq. (2.23) and Zn 
in Eq. (2.26). The expression for these operator renormalization constants 
is {Z — Zq or Z — Zii) 



Z^l + a,S, 



£ 



Z2_ 

2e 



(4.5) 
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Here Zi = 2 f3o for both operators. Further we have Z2 = 4:(3i for 0{y) in Eq. 
(2.23) and Z2 = 2 Pi for Oi{y) in Eq. (2.26). If we express this in the bare 
couphng constant it becomes 



e 



+ ciiS 



1 /I 2 n \ ^2 



(4.6) 



In the case of the electromagnetic vertex the current is conserved and we 
find Zi = 0. Finally Ku in Eq. (4.1) collects all collinear and infrared terms. 
Because of the one to one correspondence between the ln(— g^/ /i^) terms and 
the pole terms 1/e it reads up to order 



2e '-^'^ 



1 



1 

As 



(1) 



o AO) 

-PO J LI 



+ 



(4.7) 



where ff'' is a vertex dependent constant. However fl^2 is vertex independent 
as we shall see later on. From Eq. (4.1) one can derive the formula for the 
logarithm of the vertex correction. First we replace the coupling constant 
«5(/^i) by asin'^) 



ttsinl) = as(/x^) ^1 - Qsin^) (3q In j 



(4.8) 



then we integrate over Subsequently we integrate over ln(— g^///^) and 
finally we replace the renormalized by the bare coupling constant 



(4.9) 



We can now integrate Eq. (4.1) over ln(— g^//x^). However in the case of non 
conserved operators one has to be careful in the determination of the pole 
part of the vertex function. First we separate the ln(— 5^///^) terms which 

are due to the non conserved operators from the rest. Next we observe a one 
to one correspondence between the pole part and the remaining ln(— g^//x^) 
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piece. Hence the general expression for the unrenormahzed vertex correction 
becomes 



2\<^/2 



£ 9 



"I : I" 



e 

.2\ I 



-q 



1 o (0) 1 /I (1) , ,3 (0)\ 1 



(0) 
i,2 



+ 77 Z2 In 



^Inl"^ 



^2' 



2 ( ^ 



1 

77 /3o ^1 In 
2 \ /X 



^,A7* 



(4.10) 



Note that are imphcit in ^■^ and they do not contribute to the ln(— 
terms. They contribute to the pole terms (l/^)'^ only. Another remark con- 
cerns the term —2 (5q ff'} in the penultimate line of the above equation. It 
vanishes in the single pole term after coupling constant rcnormalization. The 
anomalous dimensions 7)^ jj are up to a colour factor Cj vertex independent. 
Up to two-loops they are given by 



.(0) 



^67 



« = (1,9: 



Cq — Cf , 



Cg — Ca- 



(4.11) 



However, and this is very important, ^^^o vertex independent up to a 

colour factor Cj and it is given by 



Ji,2 



a 



- Po C(2) + Ca{^- 14 C(3)) +nfTf{-^^^ 



27 



.(4.12) 
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The whole vertex dependence can be attributed to the the anomalous dimen- 
sions 7^^'' and the non-pole terms f-j^. The anomalous dimensions up to two 
loops are given by 



7, 



7, 



(0) 
99 



(4.13) 



i'J = C^(3-24C(2) + 48C(3))+C^C^(^ + ^C(2)-24C(3)) 



^(1) 

'99 



„, 32 



Cl(^ + 24 



— UfTfCA-S Uf Tf Cf . 



(4.14) 



The non-pole terms ff^^ (i — q,g) can be extracted from the actual cal- 
culations. For that purpose we also include the quark vertex correction in 
the case of electromagnetism [25], [26], [27]. Starting with the latter, this 
expression equals 

lnF^.(g2) = a,S, ( ^] ' C, 



_8+6_8 + C(2) 



.2\ £ 



44 



-30C(3)-yC'(2)[ + C^Cf 



C^|(^-12C(2) + 24C(3))i-i + 29C(2) 
f 44 / 332 ^ ^ , , N 1 



+ 



+ 



4129 11 



54 
467 



C(2)-26C(3))^ 



89173 119 



648 



9 



C(2) 



44 



9 C{S) + fC{2)\+nfTfCF 



16 112 1 



( 706 4 1 7541 28 



-fc(3)- 



(4.15) 
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The non-pole terms f^^l are 



f(0) _ 



AO) 

Jq,2 



Ci.(C(2)-8), 

C| ( - - + 29 C(2) - - C'(2) - 30 C(3)) + ( - 



89173 
648 



g C(2) + yC^(2) + — C(3)J+n; 
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-fc(3)) 



(4.16) 



Next we turn our attention to Higgs boson production. Taking the logarithm 
of the expression for the scalar Higgs boson Eqs. (3.3) we obtain 



InF^^(g') = InZo + asS, 



2\ £/2 



Ca 



+ C(2) 



_ncp)-2C(3))i + ^ + fc(2).Hc(3,} 



-fc(2)-^C(3)U.,T,cJl-f + 16C(3)' 



(4.17) 



with the non-pole constants f^^l 
C = CaC(2), 

.(0) ^ ^ / 1616 10 14. 



81 



9 



+ nfTfCF{- 



67 
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+16C(3))+C1(^ + 5IC(2) + HC(3)). 



(4.18) 



Doing the same for the pseudo-scalar Higgs boson vertex (see Eq. (3.6)) we 
get 



In FA(g^) = In Zii + tts S, 



2\ e/2 



3 + 4 + C(2) 



-llC(2)-2C(3))i + II^ + f C(2) + ^C(3) 



406 



^ ^ { 16 40 /212 , ,N 1 3182 
+-/^/^^^-3i^+9i^ + (^ + 4C(2))- 



81 



-f C(2)-^C(3)|+n,r,C^ 



(0) 



142 



+ 16C(3) 



,(4.19) 



with the non-pole constants /^^ 
/<»,> = C44 + C(2)), 

4S = %T,c4-H?H-^C(2)-l|?C(3))+n,T,C.(-l^H 



+16C(3))+Cl(I!f5 + fc(2) + Hc(3)). 



(4.20) 



Notice that there isCp ^ Ca symmetry between the electromagnetic (quark) 
and the (pseudo-) scalar Higgs (gluon) vertex correction for the '^^^^^ and //_^\ 
Prom the general form of the form factor [34] and [35] one can predict the 
coefficients of the and l/e^ terms in Eq. (4.10). These can be 

inferred from the anomalous dimensions ""yK^ii Eq. (4.2) and Eq. (4.3). 
However, and this is new, the coefficient of the single pole term 1/e can 
now also be predicted. In the renormalized vertex correction it amounts to 
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1/2 7ii^'' + fL2- The form for the electromagnetic vertex function was also 
studied in [37] and [38]. The form for the quark in [37] agrees with ours 
up to the double pole term. The same is true for the form in [38], which 
is exactly equal to ours because it is derived from [35]. However the Higgs 
vertex functions in Eqs. (4.17) and (4.19) were not discussed. Also the 
structure of the single pole term of the scattering amplitude was not given 
in [37] and [38]. The latter can be found for instance in [39] and [40] but 
the relation between the quark and the gluon was not observed for the single 
pole term. Using our approach for Eq. (3.9) in [39] and for Eq. (A. 29) in 
[40] we can express their if}^^-functions in the following way 

hP = -^#-^4^ + ^r^+^C^A)C(2), 

i = q,g, C, = Cf, Cg = CA. (4.21) 

See also Eqs. (3.10), (3.11) in [39] and Eqs. (A.30), (A.33) in [40]. It is clear 
that the structure of the single pole term is now solved. With this knowledge 
and the three-loop anomalous dimension [36] we can at least predict the 
three-loop vertex function and three-loop scattering amplitudes up to the 
double pole term. 

Summarising the above we have calculated the two-loop vertex functions 
for the scalar and pseudo-scalar Higgs boson including finite terms. We found 
agreement with [10] for scalar Higgs production for Ca = 3 and Cp = 4/3 
but the pseudo-scalar Higgs vertex correction is published for the first time. 
Further we confirmed the predictions for the form factor up to second or- 
der double pole term made in [34] and [35]. However we could also extend 
them to the coefficient of the second order single pole term which equals 
1/2'^^^^ + fo^ the quark {i = q) and the gluon {i = g). 

Acknowledgement: J. Smith would like to thank George Sterman for a 
discussion of the above results. 
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